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A New Phase at Finite Quark Density from AdS/CFT 
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We explore phases of an A/" = 2 super Yang-Mills theory at finite quark density by introducing 
quark chemical potential in a D3-D7 setup. We formulate the thermodynamics of brane embeddings 
and we find that the density versus chemical potential equation of state has rich structure. This 
yields two distinct first order phase transitions in a small window of quark density. In order words, 
there is a new first order phase transition in the region of deconfined quarks. In this new phase, the 
chemical potential is a decreasing function of the density. We suggest that this might be relevant 
to the difference in sQGP-wQGP phases of QCD. 

PACS numbers: Valid PACS appear here 



Introduction: There has been much hope that one 
might be able to use AdS/CFT [1] to describe the real 
systems after certain amount of deformations. For exam- 
ple, it has been suggested that the fireball in Relativistic 
Heavy Ion Collision (RHIC) viewed as a strongly inter- 
acting system [2, 3], has been studied using dual gravity 
models [4-8]. There have been many attempts to con- 
struct models phenomenologically closer to QCD [9]. 

More recently, there has been renewed interest in A/" = 
2 super Yang-Mills (SYM) systems with quenched fun- 
damental quark flavors studied by using a holographic 
description with probe D7-branes in the AdS^ black hole 
background [10-15]. The key observation is that we have 
confinement of quarks even in the absence of gluon con- 
finement or area law [10]. The phases of this theory are 
characterized by the brane embeddings: whether the D7- 
brane touches the black hole horizon (black hole embed- 
ding) or not (Minkowski embedding). Different types of 
embedding lead to different meson spectra. 

In this letter, we explore the phases of this theory at 
finite quark density by introducing quark chemical po- 
tential along the lines of [16, 17]. We will first establish 
a clear formulation of the thermodynamics of brane em- 
beddings. We find that we need to renormalize the finite 
chemical potential due to the divergence of the thermo- 
dynamic potentials. We will also find that apart from the 
type of first order phase transition described in [10, 11] 
at zero chemical potential, there is another class of first 
order phase transition within the black hole embedding 
category: it is indicated in Fig. 1 as a hopping between 
two black hole embeddings. 

Since black hole embeddings correspond to the decon- 
fined phase, we cautiously suggest that this new type of 
first order phase transitions might be relevant to the dif- 
ference between sQGP-wQGP in RHIC experiments. In 
particular, we find that the chemical potential in this new 
phase is a decreasing function of the density. 

We emphasize that depending on whether we control 
the system by the chemical potential (grand canonical en- 
semble) or by the density (canonical ensemble), the phase 
diagram is different. In this letter, we present the analy- 
sis and the results for the system based on the canonical 



FIG. 1: Brane embedding and phase transitions. Solid lines: 
Minkowski embedding (red) to black hole embedding (blue). 
Dashed lines: Hopping from a black hole embedding to an- 
other one. 



ensemble. The details including the description of the 
other process based on the grand canonical ensemble will 
be reported in detailed publication [18]. 

The Euclidean AdS black hole metric is given by 



U 



ds^ = {f{U)dt^ + dx^) + i?2 



f{U)U' 



+ dni\,{i) 



where f{U) = 1— (Uq/U)'^. The Hawking temperature of 
this geometry is given by T = Uq/ttB? = C/o/(V2A7ra') 
where A = gyM-^c- We introduce a dimensionless coor- 
dinate ^ defined by d^VC^ = dU^/{U^f) , so that the 
bulk geometry is 



i?2 



{fiOdf + dx^ 



R 



^2 dsl, with (2) 



dsl = d^^ + ^^dflj = dp^ + p^dill + dy^ + y^dip^, 

where we have defined = -I- and p is the radius of 
the 3-sphere. The black hole horizon is located at = 1. 
The induced metric on the D7-brane is 

d'D7 = ^ {fdt' + dx') + ^ ((1 + y'^)dp' + p'dnl) ,(3) 

where y' = dpy{p). It is interesting to notice that the 
bulk metric (1) and the induced metric (3) have the same 
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Hawking temperature. This means that the bulk and the 
brane are in equihbrium. The Euchdean DBI action of 
the D7-brane in the presence of the gauge field strength 
Fpt is 

S = Nfij,r j dt(fxdpdn3^/det{G + 2na'F)=f3V3 j dpC,{4:) 
with 



Trp"iv+ J—{l+y'^)-{F/mTV, (5) 
V ^+ 

where T7 = NfN^T'^gl,j^j/32, w±(0 = li^-" and tot = 
^VXT (the Chern-Simons term vanishes in the present 
case). 

Conserved charge and equations of motion: 

Since £ does not depend on Aq explicitly, its conjugate 
momentum is a conserved quantity: 



dFpt 

in terms of which we can write 



-<3, 



(6) 



Fpt = totQw_V(i + + (7) 

where Q = Since we have a constraint (6), to 

obtain a Lagrangian for y, we should not substitute Fpt 
into the original Lagrangian (5). The correct procedure 
is to perform a Legendre transformation 



(8) 



and then impose the conservation equation (6) to elimi- 
nate the electric field completely. The resulting "Hamil- 
tonian" is given by 



H = T{y,p)^/l^, T{y,p) = rrJ^{Q^+u,%p^).{9) 

V ^+ 

We can take this Hamiltonian as our effective Lagrangian 
for y and T may be regarded as the effective tension of 
the D7-brane. Resulting equation of motion is 



y" 



l + y'- 



d\ogT , d\ogT 
-y — = 0. 



dp 



dy 



(10) 



One can check that both the original and the effective 
Lagrangian give the same equation of motion for y. 

Chemical potential in gravity dual: In Refs. 
[16, 17], the quark chemical potential was introduced as 
the value of Aq on the D7-brane worldvolume. Here we 
define the chemical potential in a gauge-invariant fashion: 



P : 



dpK 



pt 



lim Aq. 



(11) 



For the last equality, we need to gauge fix = and 
set Ao\p^^^ = 0, which agrees with [16, 17]. 



Notice that p = J Fpt is the work to bring a unit charge 
from the UV region {p = 00) to the IR region against 
electric field Fpi. This definition (11) agrees with our 
intuition of the chemical potential as work done to add 
a quark to the system. 

For the Minkowski embedding, pmin — 0. For the black 
hole embedding, pmin = cos 9 where 9 is the angular 
coordinate on the y-p plane. 

Thermodynamic potentials: A generic grand 
potential (density) is defined by e~'^^=*^('^^ = 
j^g-pVsiH-nN) _ Here we identify the DBI action 
which is a functional of ^0 as the grand potential 
Sdbi = Then, integrating the Legendre 

transformation (8), 

j dpH = j dpC - j dpHAo^o- (12) 

Using the fact that IIaq = — Q is a constant (in p evolu- 
tion), one can rewrite the above as 



F{Q) = n{p) + pQ. 



(13) 



It is remarkable that the Legendre transformation in the 
bulk classical field theory is reinterpreted as the Legen- 
dre transformation between the canonical and the grand 
canonical ensembles in the boundary thermodynamics. 

The chemical potential enters the Hamiltonian density 
of the gauge theory at the boundary as a coupling to the 
baryon number density: 



(14) 



Therefore Q which has been originally defined as a first 
integral of the DBI action should be identified as the 
number density of quarks/baryons. Notice that the ef- 
fective tension of the brane increases as we increase the 
quark density which may have been expected. More pre- 
cisely, after considering various scale factors, we have 



(15) 



Renormalization of chemical potential: The 

Helmholtz free energy F and the grand potential are in 
fact divergent quantities since they contain a divergent 
p integral. Therefore we need to regularize them. We 
choose to subtract Fq{Q), the value of F for the D7- 
brane configuration that touches the black hole on the 
equatorial plane (y = 0). This is like a Pauli-Villars reg- 
ularization in the brane setup. We call this brane as the 
reference brane. So the renormalized free energy which 
we will calculate is defined by 



Fr{Q)^F{Q)-Fo{Q). 



(16) 



This has a far reaching effect to the chemical potential. 
To see this notice 

FR{Q) = n{p)-n{po) + {p-po)Q:=^lR{pR)+pRQ,{n) 
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where we have used the fact that Fo{Q) = ri(^o) + fJ-oQ- 
Notice that /ip is the chemical potential at the reference- 
brane configuration. The point is that when we deal with 
the renormalized quantities like Fji and Qr, we also have 
to use the renormalized chemical potential, although fi it- 
self is finite quantity from the beginning. In all numerical 
analysis, we need to use fij^ for the chemical potential. 

From now on, we delete subindcx R unless it is confus- 
ing. 

The phase structure of D3-D7: In this article, we 
study the system based on the canonical ensemble, where 
the number density Q is continuous while the chemical 
potential is allowed to jump. 

Let us expand y and Aq in the form oi y{p) — L + + 

0{p-^), Ao{p) = f.+ ^+Oip-"). We have L = 

c — c 'Z^^™^ where niq is a quark mass and c — {"ipip) 
[14, 15]. One can interpret L as the quark mass at fixed 
temperature or as inverse temperature at fixed quark 
mass. We also find that c' = — i^^xQ xhe standard 
AdS/CFT dictionary establishes that we have two pairs 
of conjugate variables, namely (c, L) and {Q,fj,). 

To analyse the phase transition, we plot the free en- 
ergy F/T"^ as a function of L for a given Q (Fig. 2) and 
track the least free energy configuration. In Figure 2, 
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FIG. 2: Free energy vs. L for fixed Q. The second phase 
transition is indicated by small kink near U'. 

we start from P' on the Minkowski branch. As L de- 
creases, the free energy decreases until we intersect the 
black hole branch at P. For smaller values of L, the black 
hole configurations have lower free energy. At P there- 
fore, the D7-brane jumps from a Minkowski embedding 
to touch the black hole horizon. This is the same type 
of first order phase transition first found in [10] for zero 
chemical potential. Decreasing L further, decreases the 
latitude (6) of the intersection of the brane with the hori- 
zon smoothly. For large Q, the story ends here. 

For very small Q, however, as we decrease L further, 
initially the latitude of intersection of brane-black hole 
goes down smoothly. But, at a critical latitude 9i, the 
embedding suddenly jumps to a smaller latitude 6*2 as 
indicated schematically in Fig. 1 (the kink near U' in 
Fig. 2). An embedding with touching latitudes between 
these two values has higher free energy as seen in Fig. 3, 
and hence is never realized. This phenomenon happens 
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FIG. 3: Free energy vs. L for fixed Q around the second phase 
transition point. 



only at finite density within a small density window, 

log Qi, -6.812 < log Q < -4.726 = log Q2* • (18) 

The full phase diagram therefore looks as in Fig. 4. No- 
tice that the horizontal axis is logQ. The chemical po- 
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FIG. 4: Phase diagram in T-Q plane. 



tential j^umps across the phase boundary. 

For Q > Q2*, there is no kink and hence the second 
phase transition disappears as we can see in Fig. 5. Fur- 
thermore, we have a second order phase transition at 
Q = Q2*- The slopes of the phase boundary hues are 
discontinuous at Q = Qi,, where three lines meet. 

T 

0.77r 
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FIG. 5: T-Q phase diagram around the new phase boundary. 

Equations of state: The relations between the 
variables (c, L) (equivalently (c, m,)) and (Q,/x) can be 
considered as the diagram of equation of state, which is 
much like the P-V diagram of Van der Waals in liquid- 
gas phase transition. We may determine these by using 
the thermodynamic relations 

_ dF OF 
dmq' ^ dQ 
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Numerical analysis shows that the relationship be- 
tween /i and Q given by Q = — is not monotonia, 
further /j, is not a single- valued function of Q. 

In Fig. 6 and 7, we show the relation between c and L 
for a representative value of Q and the relation between 
Q and /i for a representative value of L. The equations 
of state represented in these diagrams show that we have 
much richer structure than Van der Waals P-V diagram. 
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FIG. 6: /i-Q relation: We plot ^i/T vs. Q/T^ 



Discussion: Although characterized by nonstandard 
behaviour of the chemical potential, the nature of the sec- 
ond phase is not very clear from the gauge theory point 
of view. However, since both phases belong to the black 
hole embedding corresponding to deconfined quarks, it 
might be relevant to the famous difference between sQGP 
and wQGP discussed within RHIC physics [2]. In fact 
if we take the temperature and Q of the relevant re- 
gion to be 200MeV and 5 x 10^'^ respectively and choose 
iV/iVcVA/16 0(1), then the density Q - Ifm"^, which 
is in the RHIC ballpark. 
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FIG. 7: c-L relation: Here we interpret L as quark mass. 



We track the phase diagram from M' in Fig. 6 at a 
fixed value of L which may be thought of inverse temper- 
ature (or quark mass). The first phase transition takes 
place at a certain critical value of Q at point P'. Here, 
the brane embedding jumps from a Minkowski embed- 
ding to a black hole embedding which results in a jump 
of the chemical potential from P' to R' on the diagram. 
If we increase Q further, we have a second phase transi- 
tion which is realized as a jump from S' to U'. We have 
also indicated the Maxwell's construction, which allows 
us to determine the location of the phase transition (i.e., 
equality of the area of the shaded region on either side 
of a particular jump). For larger values of L however, we 
have only a single phase transition. 

The point to be noted is that in a small window be- 
tween R' and S', the chemical potential is a decreasing 
function of the density (in the deconfined phase). One 
way to achieve this is to have bound states. If this is 
true, our new phase is likely to have many bounded pairs 
of quarks. 

Referring to the c-L equation of state in Fig. 7, the two 
phase transitions are as labelled. In this case however, 
the chiral condensate decreases uniformly as a function 
of L (excepting of course, for the jumps). For large quark 
mass L, it seems to vanish as expected. 



An immediate question is the universality of this sec- 
ond phase. Since most of the properties of the branes, 
for small Q, are determined in the neighborhood of the 
horizon, we might expect that this phenomenon persists 
in (perhaps more realistic) models so long as we have a 
black hole. 

One can also look at a process which corresponds to a 
horizontal Maxwell's construction in ^-Q diagram. This 
is a process where the chemical potential is used as the 
control parameter. Such a process is more typical in the 
literature [19]. Some other related questions are about 
the effect of the density and temperature on the meson 
spectrum and on the heavy quark potentials. A limita- 
tion of this model is the absence of gluon confinement. 
We can also ask whether what we found in this paper is 
universal feature of black hole embeddings. These issues 
are currently under investigation [18]. 
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